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LONG-TERM  GOALS 

The  long  term  goal  of  this  project  is  to  determine  optimal  sampling  strategies  for  drifting  buoys,  in 
order  to  enhance  prediction  of  particle  motion  in  the  ocean,  with  potential  applications  to  ecological, 
search  and  rescue,  floating  mine  problems,  and  design  of  observing  systems. 

OBJECTIVES 

To  develop  methods  for  the  objectives  of  ODD  AS  DRI,  which  is  to  preserve  the  sampling  volume  of 
an  acoustic  sensor  array,  and  replenish  when  gaps  in  the  observational  volume  occur  due  to  the 
dispersion  of  the  sensors  with  the  underlying  Eulerian  flow  field. 

APPROACH 

The  work  is  based  primarily  on  stochastic  and  dynamical  models  of  particle  motion  and  data 
assimilation  strategies.  It  also  involves  the  use  of  coastal  circulation  models  and  processing  of 
oceanic  data  such  as  drifter  positions,  ocean  surface  currents,  and  wind  field. 

WORK  COMPLETED 

1)  Acceptance  of  a  paper  on  parameterization  of  submeso-scale  eddy  fields  as  observed  by 
high-frequency  coastal  radars  using  a  stochastic  Eulerian  velocity  model  (Caglar  et  ah,  2006). 

2)  Acceptance  of  the  Lagrangian  Analysis  and  Predictability  of  Coastal  and  Oceanic 
Dynamics  (EAPCOD)  book,  which  was  edited  by  all  Pis  of  this  project  as  well  as  by  D. 
Kirwan  (U.  Delaware)  and  T.  Rossby  (URI),  by  the  Cambridge  University  Press.  The  book  is 
in  production  to  be  published  in  March  2007. 
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3)  Participation  in  both  the  winter  and  summer  cruises  of  the  DART  (Dynamies  of  the 
Adriatic  in  Real  Time)  observational  program,  in  order  to  provide  guidelines  for  targeted 
release  of  drifters. 

4)  Development  and  implementation  of  a  eomputational  scheme  in  order  to  ineorporate 

the  effeet  of  unresolved  turbulent  motion  on  Lagrangian  particle  transport  (Haza  et  ah,  2006). 
Related  work  is  being  eondueted  in  order  to  quantify  the  exaet  influenee  of  turbulent  flows  on 
the  0(1  km)  on  partiele  motion  using  NCOM. 

5)  Continuation  of  work  on  stoehastie  boundary  eonditions  (Chin  et  ah,  2006). 

RESULTS: 

1)  DART  and  targeted  drifter  launch  strategies: 

Pis  (Griffa,  Ozgokmen)  and  postdoe  (Angelique  Haza)  have  partieipated  in  the  drifter  release 
component  of  Dynamics  of  the  Adriatie  in  Real  Time  (DART)  experiment  eondueted  jointly  by 
NATO-NURC  (led  by  Miehel  Rixen)  and  NRL  (led  by  Jeff  Book).  DART  was  an  extensive 
observational  and  modeling  experiment.  Our  partieipation  eoneemed  mainly  the  determination  of 
drifter  release  loeations  for  best  sampling  of  the  dynamieal  features.  The  primary  methods  to 
determine  the  launeh  loeations  relied  on  2-day  hindcasts  and  2-day  foreeasts  from  NCOM,  the  outputs 
of  whieh  are  provided  to  us  by  Paul  Martin  from  NRL.  This  model  has  very  high  resolution  (1  km,  1 
hour)  and  is  driven  by  4-km  resolution  COAMPS  atmospherie  foreing,  tides  and  rivers.  The  main  goal 
of  our  drifter  launches  during  the  winter  (Mareh  2006)  eruise  has  been  to  pinpoint  the  loeation  of  a 
prominent  hyperbolie  point,  which  arises  due  to  vortex  stretching  over  the  Palagruza  Sill  in  the  middle 
of  the  Adriatie.  Sueh  hyperbolie  points  ean  faeilitate  high  levels  of  partiele  dispersion,  and  thus  are 
useful  to  inerease  the  sampling  coverage.  While  the  eoneept  of  hyperbolic  points  as  prominent 
features  eontrolling  Lagrangian  dispersion  is  well  developed  in  theory  and  numerieal  models,  their 
existenee  and  influenee  in  reality  requires  more  support  from  observations.  During  the  seeond  cruise 
(August  2006),  the  main  effort  was  spent  on  eapturing  the  reeireulating  eddies  forming  downstream  of 
the  Gargano  Cape.  Sueh  eddies  appear  to  form  mostly  in  the  summer  season  and  have  a  signifieant 
impaet  on  eoastal  regions  by  entraining  water  masses  off  the  eoastal  shelf.  Two  methods  have  been 
used.  First,  historieal  data  sets  have  been  analyzed  based  on  drifters  released  by  P.  Poulain  in  the 
Adriatie  Sea  over  the  last  deeade.  Seeond,  the  interseetion  points  of  finite-seale  Lyapunov  exponents 
(FSLEs)  from  2-day  baekward  and  2-day  forward  particle  integrations  have  been  traced.  We 
emphasize  that  FSLEs  are  not  only  useful  to  deteet  hyperbolie  points  and  shear  zones,  which  are  useful 
to  enhanee  eoverage  and  independent  sampling,  but  also  zero-ESEE  regions  are  suitable  for  transport 
of  partiele  elusters  without  deformation  and  dispersion,  whieh  is  one  of  the  optimization  eriteria  of 
ODD  AS.  To  our  knowledge,  this  is  the  first  time  that  targeted  drifter  launch  strategies  have  been 
employed  real-time  during  an  observational  program. 

We  have  eondueted  a  preliminary  analysis  of  data  eolleeted  in  the  winter  eruise  in  eomparison  to 
predietions  based  on  NCOM  fields.  During  this  experiment  (Mareh  11-23,  2006),  a  total  of  17  drifters 
were  launehed,  of  whieh  12  have  survived  for  a  period  longer  than  2  days.  Trajectories  of  all  these 
drifters  are  shown  in  Eig.  1.  They  provide  a  good 
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Figure  1:  Trajectories  of  DART 
drifters  during  March  11-30. 


X  (16  ±  1,  r=5),  launched  Mar  16, 17:00 


X  (19  ±  1,  r=5),  launched  Mar  19,  20:00 
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Figure  2:  Comparison  of  the  real  DART  (green  and  pink)  and  synthetic  NCOM  (black)  trajectories 
from  the  March  2006  experiment.  Two  drifters  were  launched  based  on  finite-scale  Lyapunov 
exponents  (plotted  in  the  background)  from  the  NCOM  daily  output,  and  constraints  due  to  the 
predetermined  cruise  path  (regular  gray  line  in  the  background).  Note  that  clusters  of  synthetic 
trajectories  provide  a  reasonable  prediction  of  the  real  surface  drifter  motion. 
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coverage  of  the  region,  including  the  boundary  eurrent  and  the  interior,  which  indicates  that  one  of  the 
overall  goals  of  the  experiment  was  achieved.  In  order  to  examine  the  whole  data  set,  we  re-initialized 
daily  the  synthetic  trajectories  around  the  observed  drifter  positions  and  integrated  forward  for  7  days. 
Daily  initialization  is  motivated  by  the  fact  that  the  Lagrangian  time  scale  in  this  region  is 
approximately  1-1.5  days,  so  that  results  on  Lagrangian  predictability  are  expeeted  to  be  independent 
over  1  to  2  day  intervals. 

By  far  the  most  striking  result  from  this  analysis  is  that  NCOM  can  satisfactorily  describe  surface 
Lagrangian  transport  (Fig.  2)  in  the  Adriatie  Sea,  which  is  very  promising  for  operational  programs. 
The  summer  experiment  involved  rather  dramatic  changes  between  Bora  (recireulating  eddy 
suppressing)  and  Sciroeco  (recirculation  eddy  favoring)  winds.  Overall,  DART  experiment  was  an 
exeellent  opportunity  not  only  for  testing  Lagrangian  methods  in  real-time  applications,  but  also  for 
establishing  new  collaborations  with  scientists  specializing  in  operational  oeeanography  (M.  Rixen,  P. 
Martin,  J.  Book)  and  maintaining  researeh  with  others  (A.  Moleard,  P.  Poulain).  We  envision  to 
publish  two  joint  papers  from  these  experiments,  a  short  deseriptive  one  to  be  prepared  in  2006,  and  a 
comprehensive  paper  on  the  Lagrangian  eomponent  in  2007.  In  addition  to  Lagrangian  researeh,  one 
of  the  byproducts  of  our  DART  participation  is  that  other  in-situ  and  remote  sensing  observations  will 
be  useful  for  the  Ph.D.  dissertation  of  M.  Marcello  on  the  impact  of  capes  on  coastal  circulation,  and 
M.  Rixen  will  serve  in  his  eommittee. 

2)  Lagrangian  subgridscale  model  for  transport  improvement: 

An  aecurate  estimation  of  Lagrangian  transport  in  the  oeean  is  important  for  a  number  of  practieal 
problems  such  as  dispersion  of  pollutants,  biological  species,  and  sediments.  Forecasting  of  the 
Lagrangian  pathways  necessarily  relies  on  the  aecuracy  of  ocean  and  eoastal  models.  However,  these 
models  include  a  number  of  errors  due  to  forcing,  missing  physics  and  unresolved  turbulenee  that 
propagate  directly  from  the  Eulerian  velocity  field  to  the  Lagrangian  transport. 

In  this  study  (Haza  et  ah,  2006),  a  method  is  developed  in  order  to  reduee  errors  in  Lagrangian 
transport  by  employing  the  statistieal  drifter  behavior  estimated  from  observational  data  sets.  The 
method  is  shown  to  work  well  using  both  a  so-ealled  Markov  veloeity  field  model,  representing  an 
idealized  turbulent  flow  field,  and  in  the  eontext  of  NCOM  configured  in  the  Adriatie  Sea  for  realistie, 
high-resolution,  complex  oeean  flows 
(Fig.  3). 

The  simplicity  and  computational  efficiency  of  this  technique,  eombined  with  applieability  to  oeean 
models  at  a  wide  range  of  resolutions,  appears  promising  in  light  of  the  ehallenge  of  eapturing  exactly 
the  oeeanic  turbulent  fields,  whieh  is  critical  for  Lagrangian  dispersion. 

3)  Impact  of  small-scale  turbulent  flows  on  particle  motion: 

Traditionally,  the  impaet  of  small-seale  turbulent  flows,  sueh  as  inertial  oscillations, 
divergence/convergence  zones  near  coastal  flows  generated  by  upwelling/downwelling,  strong  wind 
events  and  vertical  mixing  have  been  neglected.  This  is  not  only  due  to  inadequate  spaee/time 
sampling  of  the  observations,  but  also  because  Lagrangian  methods 
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Figure  3:  Fifteen-day  long  trajectories  of  121  particles  from  Adriatic  NCOM  in 
various  experiments  where  the  parameters  of  the  Lagrangian  subgridscale  model 
are  varied  to  modify  the  model-generated  trajectories  (middle panel)  to  be  much 
less  (upper  left  panel)  or  much  more  (lower  right  panel)  energetic 

(Haza  et  ai,  2006). 
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based  on  dynamical  system  theory  often  work  well  only  in  flow  regimes,  in  which  the  Lagrangian  time 
scale  is  smaller  than  the  Eulerian  time  scale.  In  light  of  very  high-frequency  sampling  (e.g.,  GPS)  data 
sets  becoming  available  for  oceanic  community,  high-resolution  numerical  models,  highly  variable 
nature  of  coastal  surface  flows  (e.g.,  NCOM  and  DART  experiment),  the  question  arises  how  metrics 
quantifying  dispersion  of  particles,  such  as  FSLEs,  are  sensitive  to  time/space  filtering  of  scales.  If 
they  are  sensitive,  what  can  be  done  to  compensate  for  the  dynamics  filtered  out  during  such 
space/time  averaging?  These  problems  are  being  pursued  with  A.  Poje  (CUNY)  based  on  high- 
resolution  outputs  of  NCOM  in  the  Adriatic.  Our  initial  finding  confirms  that  time  smoothing  changes 
the  behavior  of  particle  dispersion  (Fig.  4).  The  impact  of  space  smoothing  will  be  tackled  next. 


Figure  4:  Dependence  of  FSLE  on  particle  separation  distance  estimated  from 
NCOM  configured  in  the  Adriatic.  This  plot  shows  that  low-pass  filtering  with  1.5 
and  3-day  intervals  (red  and  green  curves,  respectively)  leads  to  significant  errors 
in  submesoscale  relative  dispersion  behavior,  but  large-scale  dispersion  is 

minimally  effected. 


IMPACT/APPLICATIONS 

The  investigation  of  the  predictability  of  particle  motion  is  an  important  area  of  study,  with  a  number 
of  potential  practical  applications  at  very  different  scales,  including  searching  for  persons  or  valuable 
objects  lost  at  sea,  tracking  floating  mines,  ecological  problems  such  as  the  spreading  of  pollutants  or 
fish  larvae,  design  of  observing  systems  and  navigation  algorithms. 
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